Configuration interaction calculation of hyperfine and P,T-odd constants on 207 PbO 
excited states for the electron EDM experiments 
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We report first configuration interaction calculations of hyperfine constants A\\ and the effective 
electric field Wd acting on the electric dipole moment of the electron, in two excited electronic states 
of 207 PbO. A new combined scheme of correlation calculations is also developed and first applied 
to studying the core properties. The obtained hyperfine constants, Aii = —3826 MHz for the a(l) 
state and A\\ = 4887 MHz for the B(l) state, are in very good agreement with the experimental 



data, -4113 MHz and 5000 ± 200 MHz, respectively. We find W d = 
o(l), and W d = - (8.0 ± 1.6) -10 24 Hz/(e ■ cm) for B{1). 
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Introduction. The search for the electric dipole mo- 
ment (EDM) of the electron, d e , remains one of the most 
fundamental problems in physics. Up to now only up- 
per limits were obtained for \d e \; the latest constraint, 
\d e \ < 1.6 ■ 10 -27 e-cm, was obtained in an experiment 
on atomic Tl [lj. The search for d e in heavy-atom po- 
lar molecules was initiated by theoretical investigations 
0, HI on radicals with one unpaired electron. (The first 
experimental bound on d e using such a system was re- 
cently obtained, using YbF 0.) It was since noted that 
experiments on the excited a(l) [j| or 5(1) ||| states of 
PbO, which have two unpaired electrons, may make it 
possible to search for d e with sensitivity 2-4 orders of 
magnitude higher than the current limit. An important 
feature of such experiments is that the knowledge of the 
effective electric field, Wd, seen by an unpaired electron 
is required for extracting d e from the measurements. Wd 
can not be obtained in an experiment; rather, electronic 
structure calculations are required for its evaluation. Re- 
cently, a semiempirical estimate of Wd for the a(l) state 
and ab initio relativistic coupled cluster (RCC) @ 
calculations with the generalized relativistic effective 
core potential (GRECP) |lcL fll| of Wa for both the o(l) 
and B(l) states were performed. The nonvariational one- 
center restoration (NOCR) technique [Tel HI HI lli| 
was employed for obtaining proper electronic densities 
close to nuclei after GRECP /RCC calculations of heavy- 
atom molecules. The semiempirical estimate of Ref. [jj, 
\Wd\ > 12 x 10 24 Hz/(e-cm), is three times higher than 
the RCC result (see Table QJ. Calculations performed in 
demonstrated the need for a higher level of accounting 
for correlation in the valence region for the excited a(l) 
and -B(l) states. The main problem was that the Back- 
space RCC with single and double excitations (RCC-SD) 
version used in was not optimal for accounting for the 
nondynamic correlations (see 0] for details of RCC-SD 
calculations of the Pb atom) , though the potential of the 
RCC method for electronic structure calculations is very 
high in prospect The estimated error bounds put 



the actual Wd value between 75% and 150% of the value 
calculated by the GRECP/RCC/NOCR approach. 

The main goal of the present work is to develop a 
method of calculation which could provide substantially 
higher accuracy and reliability for Wd in the excited 
states of PbO and could be efficiently applied to other 
heavy-atom molecules. As an accuracy check for the cal- 
culation of the electronic structure near the Pb nucleus, 
the hyperfine constants (An ) are also calculated. For this 
purpose in the present work the SODCI code 0] for 
the spin-orbit direct configuration interaction (CI) calcu- 
lation is employed, in which the relativistic (spin-orbit) 
scheme and new criteria Ell of the configuration 
selection proposed by us earlier are incorporated. New 
codes for calculation of a one-electron density matrix 
with the CI wavefunction and restored four-component 
molecular spinors have been developed. The 10-electron 
GRECP/SODCI/NOCR calculations are performed for 
Ay and Wd in the a(l) and B(l) states of PbO with 
wavefunctions restored in the Pb core (see 0, H3, for 
recent two-step calculations). 

The space of a SODCI calculation consists of a set of 
many-electron configurations obtained by selecting the 
singly- and doubly-excited configurations with respect 
to some multiconfigurational reference states when using 
some selection thresholds (T) [Isl 122^ . Because of lim- 
ited computational resources, all the possible singly- and 
doubly-excited configurations (without selection) cannot 
be usually included to the SODCI calculation. The cor- 
rections accounting for zero threshold (T=0) [12] and 
higher-order excitations ( "full-CI" ) allow one to esti- 
mate the contribution to the total energy of all the con- 
figurations omitted in the SODCI calculation (see below 
for more details). In the present work we first apply the 
T=0-, full-CI- and "outercore correlations" -type correc- 
tions to calculating the properties which cannot be ob- 
tained from the energy curves (surfaces). The applicabil- 
ity of these corrections is justified and discussed in details 
in [24|. The developed method has allowed us to reduce 
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the deviation of the calculated An from its experimental 
value for a(l) in about five times, from 34% to 7%, as 
compared to the previous calculation |jj . That deviation 
for 5(1) is reduced from (18 ± 4)% to (2 ± 4)%, i.e., it 
is now within the experimental uncertainty ±4%. The 
SODCI results without applying the T=0, full-CI and 
outercore correlation corrections give the deviations 12% 
for a(l) and (8 ± 4)% for B{1). As a consequence, sub- 
stantially more reliable values for Wd are obtained. Our 
present error estimate on Wd is relatively more conserva- 
tive than in since \Wd\ is increased from the RCC-SD 
datum on 50% for a(l) and decreased on 18% for B(l) 
(these notable changes are, however, within the error lim- 
its declared in ||). The expression for Wd is given by 
Refs. jElSEEl 

W d = ^- e ^x\Y, H d(^x) , (1) 

i 

where ^ x is the wavefunction for a(l) or B(l) state, 
and il = (*S/x\J ■ n\^>x) = ±1, J is the total electronic 
momentum, n is the unit vector along the molecular axis 
directed from Pb to O, 

Hd = 2 de (° Q (2) 

E is the inner molecular electric field, <x are the Pauli 
matrices. The hyperfine constant A\\ is determined by 
the expression |27j 

-4,,4^<**iE(^)zi**), m 

where I and /ipb is the spin and magnetic moment of 
207 Pb, a.i are the Dirac matrices for the i-th electron, and 
Ti is its radius-vector in the coordinate system centered 
on the Pb atom. 

Methods and calculations. The above properties were 
calculated in two steps. At the first step, a 22-electron 
GRECP for Pb |2^| simulating interaction with the ex- 
plicitly excluded Is to 4/ electrons is used. In addi- 
tion, the 5s 2 5pf^ 2 5py 2 5dy 2 5d^ 2 shells of lead and the 
Is 2 shell of oxygen were frozen (see ^{j for details) and 
the residual ten electrons were treated explicitly in the 
subsequent molecular calculations. The basis set on Pb 
(15sl6pl2d9/)/[5s7p4d2/] as well as GRECP are the 
same as those used in paper |9j. The basis set was op- 
timized for calculation of properties determined mainly 
by the electronic wave function at the Pb core. The 
description of the basis set generation procedure can 
be found in Refs. (lil Ej| . The correlation-consistent 
(10s5p2dl f)/[4s3p2 dl f] basis of Dunning listed in the 
MOLCAS 4.1 library |3JJ] was used for oxygen. 

The leading AS coupling terms and configurations 
for the a(l) and B(l) states are 3 E + ofcrf cfTr 3 ^ and 



3 IIi (T 2 o'2 cr 3 7r i 7r 2! respectively. The molecular orbitals 
used in the CI calculations are obtained by the restricted 
active space self consistent field (RASSCF) method [111 

twith the spin-averaged part of the GRECP (AGREP) 
, i.e. only scalar-relativistic effects are taken into ac- 
count in the RASSCF calculation. Because the a(l) state 
is of primary interest, molecular orbitals were generated 
specifically for the lowest 3 S + state. This set of orbitals 
was used for the subsequent CI calculations of both the 
a(l) and B(l) states. In the RASSCF method, orbitals 
are divided into three active subspaces: RAS1, with a 
restricted number of holes allowed; RAS2, where all pos- 
sible occupations are included; and RAS3, with an upper 
limit on the number of electrons. In this calculation, 
no more than two holes in RAS1 and two electrons in 
RAS3 are allowed. Using the C 2v point group classifi- 
cation scheme, two Ai orbitals in RAS1, six orbitals in 
RAS2, (two each in Ai, Bi, and B 2 irreps) and 50 (20 Ai, 
6 A2, 12 Bi, and 12 B2) in RAS3 subspaces are included. 

Next the spin-orbit CI approach with the selected 
single- and double-excitations from some multiconfigu- 
rational reference states ("mains") j3^ is employed on 
the sets of different AS many-electron spin- and space- 
symmetry adapted basis functions (SAFs). In the SODCI 
code, the double C 2v group, C 2v , is used to account for 
the spin and space symmetry of the PbO molecule, in- 
stead of the more restrictive symmetry group C^ v , which 
could in principle be employed. In the C 2v classification 
scheme, the doubly-degenerate a(l) and B{\) states have 
the components in the B* and B 2 irreducible represen- 
tations (irreps). The operators Wd and An have nonzero 
matrix elements only between the wavefunction compo- 
nents from B* and B 2 irreps. So, one must know both 
these components of the a(l) and B(l) states to calcu- 
late Wd and .A 1 1 when working in the C 2v group. The 
SAFs from the 2S+1 C2„-irreps (constructed on the ba- 
sis of AGREP /RASSCF pseudoorbitals) , singlet (^1), 
triplets ( 3 Ai, 3 A 2 , 3 B 2 ) and quintets ( 5 Ai, 5 A 2l 5 Bi, 
5 B 2 ), were included in calculations of the components be- 
longing to the BI irrep and equivalent calculations were 
performed for those laying in the B 2 irrep. (Alterna- 
tively, the components belonging to the B 2 irrep can be 
constructed by acting on those from the B\ irrep by the 
operator of projection of the total electronic momentum 
on the molecular axis.) The reference space, 
consisted of No = 2517 SAFs having the largest coeffi- 
cients in the probing CI calculation. The single and dou- 
ble excitations from this reference space produce about 
A T=0 «175 000 000 SAFs, ' 2) }£Zf°- The molecular 
Hamiltonian H in our calculations is presented as 

H = H° + V , (4) 

where H° is an unperturbed spin-averaged Hamiltonian 
constructed to be diagonal in the given many-electron 
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basis set: 



where n = i,j and 
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and V is a perturbation that includes the two-electron 
operator describing correlations and a one-electron effec- 
tive spin-orbit operator. Only the most important SAFs 
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were selected (with all SAFs from the set) for 

the subsequent CI calculation. In the above equation, 
£ x and are the eigenvalue and eigenfunction of H 
for a(l) or B(l) state in the subspace About 
120000, 500 000, 1 100 000, and 2 000 000 SAFs were se- 
lected (see Table HJ for the thresholds (in 10~ 6 a.u.) 
Xi=0.1, T 2 =0.01, T 3 =0.0025 and T 4 =0.001, respectively. 
(Selection on the CI coefficients obtained with the first- 
order perturbation theory [Hj | was also used in our cal- 
culations of Ay and Wd but the obtained results for the 
same number of selected SAFs were very close to those 
obtained using © and they are not presented in the pa- 
per.) The final wavefunction is written as 



*2 = ££<? 

71=0 J 



TiW s (n) 
XJ V J ' 



(7) 



where summation is performed over selected SAFs for a 
given threshold Ti and those from the set. With 

this wavefunction, energy and any property W can be 
easily calculated: 
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Then, the linear T=0 |22j and generalized Davidson (full- 
CI or FCI) @ corrections are applied to the calculated 
properties. The T~0 correction extrapolates the results 
of calculations with two different thresholds Ti and Tj to 
a result of the calculation with all singly- and doubly- 
excited SAFs included. The full-CI correction approxi- 
mates the effect of all possible excitations of higher level 
(for a given number of correlated electrons and fixed basis 
set). The T=0 correction for energy is given by 
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The full-CI corrected energy is given by 
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where I ^ | 2 are weights of the reference SAFs in 



*jf ln , T min = min(T l ,r j ). If all singly-excited SAFs, 
{$^}, are included in calculation and the number of 
the reference SAFs is large enough, then similar expre- 
sions for T=0 and FCI corrections (in more details they 
will be discussed in |2^|) for the one-electron property W 
can be written: 
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(*^-|W|*^). Though equations JTSJ and 
(|14|l look similar to (|10fl and l|12fl . they have essentially 
different physical sense. More details about the features 
of constructing the reference space and selection proce- 
dure are given in Refs. [II HI. 

Before calculating A\\ and Wd, the shapes of the four- 
component molecular spinors are restored in the inner 
core region after the two-component GRECP calculation 
of the molecule. For this purpose the NOCR method 12j 
is applied at the second step of the calculation. (See 2lJ 
for the currently used formulation of the NOCR scheme 
and || for details of restoration of the four-component 
molecular spinors near the Pb nucleus). 

We designate the 5s, 5p, 5d orbitals of lead and Is or- 
bital of oxygen as the " outercore" , and the o~i , o~2 , 03 , 
7i"i, 7r 2 orbitals of PbO (consisting mainly of 6s, 6p or- 
bitals of Pb and 2s, 2p orbitals of O) as valence. In the 
CI calculations we take into account only the correla- 
tions between valence electrons. At the final stage of our 
calculation, we estimate the contribution from correla- 
tions between valence and outercore electrons (including 
high-order correlations between outercore electrons) as 
difference in the results of the corresponding 30- and 10- 
electron GRECP/RCC calculations. (See also [15j where 
this correction is applied to the Pb atom) . Such correla- 
tions are designated in the text and tables as "outercore 
correlations" . 

Results and discussion. Calculations were performed 
at two internuclear distances, R = 3.8 a.u. (as in the 
RCC calculations of 9]), and R — 4.0 a.u. The latter 
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point is closer to the equilibrium distances both for of 1) 
(R e = 4.06 a.u. 33]) and for B(l) (R» = 3.914 a.u. [Ij). 
However, in the RCC calculations Q, R = 3.8 a.u. was 
used because of a problem with convergence at R = 4.0 
a.u. The calculated values with the one-center expansion 
of the molecular spinors in the Pb core on either s, s;p 
or s;p; d partial waves are collected in Tabled The final 
data are obtained as the result of linear extrapolation to 
the experimental equilibrium distances. 

As was noted in |9j when comparing the RASSCF re- 
sults with those obtained by the RCC method, the spin- 
orbit interaction changes An and Wd dramatically, so 
even the sign for Wd is changed. One can see a simi- 
lar picture from comparison of the RASSCF and CI data 
obtained in the present work. The second point to notice 
is a significant difference of the results obtained at the 
internuclear distances R = 3.8 and R — 4.0 a.u., espe- 
cially, for the Wd parameter. It is increasing for the o(l) 
state and is decreasing for the B(l) state, by 15% and 
30%, respectively. The agreement with the experimental 
datum attained in the CI calculation of An for the a(l) 
state at the point R = 4.0 with the threshold T 4 =0.001 
is 13%. After applying T=0 and FCI corrections and 
taking into account the outercore correlations the agree- 
ment is improved to the level of 9%. The calculated An 
value for the B(l) state coincides with the measured da- 
tum within the experimental uncertainty of 4%. The 
results with the T=0 correction are close to those with 
the smallest threshold T\ =0.001. To check the reliabil- 
ity of the linear T—0 correction j2^| for An and Wd we 
have calculated it for three different pairs of thresholds: 
Xi and T2, T2 and T3, T3 and T4. For An all three pairs 
give the same result within 1% of accuracy. As to Wd, 
the result with the T=0 correction for the first pair dif- 
fers from those for the other two pairs by 8% for a(l) 
and by 5% for B(l). However, the last two pairs again 
give the same result within the accuracy of 1% for both 
states. So, we suggest, that the T=0 limit is determined 
for our main configurations with an accuracy of 1% for 
both considered properties. Because the reference space 
is large enough, we have a small FCI correction. When 
taking into account outercore contributions at the point 
R = 4.0 a.u. we used the results of the RCC calculation 
at the point R = 3.8 a.u. This assumption seems reason- 
able for several reasons. First, the core should relax less 
than the valence region when R is changed. In addition, 
because of the spatial separation between the core and 
valence electrons, core-valence correlation contributions 
should be more stable than valence- valence ones. Finally, 
since these contributions arc relatively small, we expect 
errors due to this approximation not to be severe. 

We next discuss the uncertainty in the calculated Wd 
values. (We confine this detailed discussion to the a(l) 
state; similar considerations are applied to B(l).) Since 
Wd is sensitive to the wavefunction and its derivative at 
the Pb nucleus, it is natural to use the value of Ay (which 



is also singular close to the Pb nucleus) as a benchmark 
for accuracy and stability. Thus, the 10% deviation of 
the calculated value of An from the experimental value 
represents an obvious lower bound for the accuracy of 
Wd- It appears that Wd is less computationally stable 
than An, however, as may be noted from the variation 
in values at various stages of the calculation, shown in 
Table 0] Thus, this simple argument may underestimate 
the error in Wd- Because of the good convergence of the 
CI calculation, it appears that the deviation in Ay can 
be due in part to sensitivity of the results to the value 
of R, and in part to incomplete accounting for outercore- 
valence correlations. Consideration of the former effect 
leads us to a possible error in Wd of ~ 10%. The latter is 
more difficult to account for numerically, but (based on 
experience from similar calculations on HgF, PbF, YbF, 
and BaF [3 HJ |H Ull ) we believe that the apparent 
insensitivity of Ay to outercore- valence correlations (com- 
pared to the case for Wd) is an artifact unlikely to per- 
sist in more advanced calculations. Finally, we note that 
every improvement to the calculation (increasing corre- 
lation threshold; inclusion of core- valence correlation; in- 
creasing value of R towards the equilibrium value) actu- 
ally increases \Wd\- Thus, the true value of \Wd\ is not 
likely to be much lower than the present calculated value. 
The linear extrapolation of the calculated A|| and Wd to 
the experimental equilibrium distance for a(l), R e = 4.06 
a.u. H3, gives us -3826 MHz and -6.1-10 24 Hz/(e • cm), 
respectively, and the deviation of the extrapolated Ay 
value from the experiment now is only 7%. Account- 
ing for the above arguments, the estimated error bounds 
put the actual Wd value between 90% and 130% of our 
final value. We obtain Wd for B(l) by linear interpola- 
tion of the data obtained at R = 3.8 and 4.0 a.u. to the 
experimental equilibrium distance, R e = 3.914 a.u. |34j |. 
yielding Wd = — 8.0-10 24 Hz/(e • cm), with an estimated 
error range of ±20%. Similar interpolation for Am on 
B(l) gives 4887 MHz, which is within the uncertainty of 
the current experimental value (3(|, 5000 ± 200MHz. 

Conclusion. A new combined method for accurate 
calculating core properties in excited states of heavy- 
atom molecules is developed. It has allowed us to ob- 
tain accurate Ay and Wd values for the a(l) and B(l) 
states of PbO in spite of their complexity. The developed 
method and codes open an opportunity of reliable theo- 
retical study of whole multitude of core properties in such 
computationally challenging objects as low-lying states of 
heavy-atom systems with valence p-electrons and dense 
spectra. 
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TABLE I: Calculated parameters A\\ (in MHz) and W d (in 10 24 Hz/(e ■ cm)) for the a(l)and B(l) states of 207 PbO at 
internuclear distances 3.8 and 4.0 a.u. The experimental values of Ay are —4113 MHz in a(l) |37l |. and 5000 ±200 MHz in 5(1) 
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Internuclear distance R 


= 3.8 
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lOe-RASSCF 


-894 


-1505 


-1503 


0.73 


0.70 








0.0 


0.0 


lOe-RCC-SD [9] 






-2635 


-2.93 
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3878 


-11.10 


-10.10 


30e-RCC-SD [9] 






-2698 




-4.10 






4081 


-9.10 


-9.70 


outer-core (30e-RCC-SD - lOe-RCC-SD) 






-63 




-1.05 






203 




0.40 


lOe-CI (reference)" 


-406 


-1877 


-1874 


-0.74 


-0.83 


731 


3785 


3805 


-7.67 


-7.17 


lOe-CI (T=0.1) 


-472 


-2930 


-2926 


-2.12 


-2.21 


393 


4051 


4074 


-9.85 


-9.40 


IOp-CT (T— 01 1 


-430 


-3222 


-3218 
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-3.13 
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(10e-CI + T=0 + FCI + outercore) 






-3509 




•5.18 






4785 


_ 


10.24 


(T is m 10 a.u.) 








Internuclear distance R 


= 4.0 


a.u. 






lOe-RASSCF 


-770 


-1384 


-1383 


1.05 


1.00 








0.0 


0.0 


lOe-CI (reference) 


-459 


-2026 


-2025 


-0.64 


-0.72 


966 


4127 


4150 


-6.69 


-6.22 


lOe-CI (T=0.1) 


-479 


-3125 


-3124 


-2.34 


-2.44 


525 


4332 


4357 


-7.79 


-7.35 


lOe-CI (T=0.01) 


-449 


-3458 


-3458 


-3.50 


-3.61 


495 


4565 


4590 


-7.38 


-6.94 


lOe-CI (T=0.0025) 


-426 


-3536 


-3536 


-3.97 


-4.08 


481 


4636 


4662 


-7.45 


-7.02 


lOe-CI (T=0.001) 


-422 


-3571 


-3571 


-4.19 


-4.31 


480 


4666 


4692 


-7.49 


-7.07 


lOe-CI + T=0 






-3625 




-4.65 






4739 




-7.15 


lOe-CI + T=0 + FCI 






-3689 




-4. 81 






4762 




-7.18 
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(lOe-CI + T=0 + FCI + outercore) 6 






-3752 




■5.86 






4965 




-6.78 



a "Reference" means that the CI calculation was performed with the reference (main) SAFs only. 

b It is assumed that the outercore contribution at the internuclear distance R — 4.0 a.u. is approximately the same as is at 
the point R = 3.8 a.u. 
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